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Abstract. We compare the photometric accuracy of spectra and images in flares 
observed with the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) space- 
craft. We test the accuracy of the photometry by comparing the photon fluxes 
obtained in different energy ranges from the spectral-fitting software SPEX with 
those fluxes contained in the images reconstructed with the Clean, MEM, MEM- 
Vis, Pixon, and Forward-fit algorithms. We quantify also the background fluxes, 
the fidelity of source geometries, and spatial spectra reconstructed with the five 
image reconstruction algorithms. We investigate the effects of grid selection, pixel 
size, field-of-view, and time intervals on the quality of image reconstruction. The 
detailed parameters and statistics are provided in an accompanying CD-ROM and 
web page. We find that Forward-fit, Pixon, and Clean have a robust convergence 
behavior and a photometric accuracy in the order of a few percents, while MEM 
does not converge optimally for large degrees of freedom (for large field-of-views 
and/or small pixel sizes), and MFM-Vis suffers in the case of time- variable sources. 
This comparative study documents the current status of the RHESSI spectral and 
imaging software, one year after launch. 
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1. INTRODUCTION 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 
was launched on 2002 Feb 5 and operated successfully through the 
first year, recording a total of ~ 7500 solar flares. At this time, after 
one year of the mission, most of the software has matured to be both 
reliable and accurate, and instrumental parameters have been improved 
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from in-flight data to great precision. Thus, it is timely to conduct a 
quantitative comparison of the performance of the various available 
imaging algorithms, to compare the spatial and spectral features in the 
reconstructed images, and to cross-calibrate the photon fluxes extracted 
from images against those from the (non-imaging) spectral fitting soft- 
ware. Such a quantitative comparison will also be useful to provide 
estimates of uncertainties and systematic errors in spectral fitting and 
spectral inversions, as well as to serve as a reference for monitoring 
future changes in the instrumental parameters and software develop- 
ments. All image reconstructions are performed here with the software 
version available in March 2003. 

The RHESSI spacecraft and its mission are described in the overview 
article of Lin et al. (2002), the RHESSI spectrometer in Smith et al. 
(2002), and the RHESSI imaging concepts in Hurford et al. (2002). 
Further technical articles that may be of interest in the context of this 
study are on RHESSI data analysis software (Schwartz et al. 2002), on 
reconstruction of RHESSI images with forward-fitting (Aschwanden et 
al. 2002a), and on pixon image reconstruction (Metcalf et al. 1996). 
We choose for this comparative study as main event the flare of 2002 
Feb 20, 11:06 UT, because this event was already subject of various 
analyses in a number of previous publications (Krucker & Lin 2002, 
Sui et al. 2002, Vilmer et al. 2002, Brown et al. 2002, Aschwanden et 
al. 2002b, Kontar et al. 2002), as well as the two additional flares of 
2002 May 18, 19:15 UT, and 2002 Jul 23, 00:30 UT. 



2. RESULTS OF DATA ANALYSIS 

The details of the data analysis are provided on the accompanying CD- 
ROM (or on the web-page http://www.lmsal.com/~aschwand/eprints/ 
2003-photo/index.htmT). Here we summarize only the main results of 
this analysis, which includes a detailed parametric study of the 2002 
Feb 20 flare (§A1-A11 on CD-ROM), as well as additional tests on flares 
with different source morphologies and flare positions, i.e., the single- 
source flare from 2002 May 18 (§B1 on CD-ROM) and the complex- 
source X-ray flare of 2002 Jul 23 (§C1 on CD-ROM). 

2.1. Spectral Photometry (§A1 on CD-ROM) 

The photometry can be obtained from two different software codes 
in the RHESSI SolarSoftWare (SSW), either from the spectral fitting 
package SPEX (Schwartz et al. 2002; Smith et al. 2002), or from the 
RHESSI image reconstruction algorithms. Both codes use the same 
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Table I. Ratios of image-integrated photon fluxes to the total spectral fluxes from SPEX. 
The last column shows the means and standard deviations of all 5 imaging methods. 
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input in form of time-tagged photon event lists, but the RESSI response 
function is applied by independent modules to obtain photon counts 
in a given energy range. A first basic photometry test that verifies the 
consistency of the software is therefore to compare the photon counts 
of the SPEX spectra in different energy ranges with those obtained 
from images in the same energy ranges. The time profile of the flare, 
the photon spectrum, the photon count rates and pre-flare background 
in 8 energy ranges are all shown in §A1 on the CD-ROM. The pre-flare 
background amounts to a few percent of the peak flux (at energies < 40 
keV). 

2.2. Image Photometry (§A2 on CD-ROM) 

We conduct a photometric test by comparing the photon rates (in units 
of photons s _1 cm' 2 ) integrated over each image with those obtained 
from the total flux spectra calculated with SPEX in identical energy 
bins. We find that these total fluxes are well conserved in most of the 
imaging reconstruction algorithms (see ratios in Table 1), they agree 
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Figure 1. Schematic cartoon of solar and non-solar sources that produce modulated 
flux F mo d(t), unmodulated flux F unm (t), and non-flare sources Fb ac k(t) in the time 
domain of the modulation time profiles. 

with SPEX within a few percents for all 5 imaging algorithms (Forward- 
fit, Pixon, Clean, MEM, MEM-Vis). The photometric accuracy seems 
to be uniform in different energy ranges (last column in Table 1). 
This test essentially verifies that the flux normalization is correctly 
and consistently implemented in the different software modules of each 
imaging algorithm. 

2.3. Background Components (§A3 on CD-ROM) 

The imaged flux consists of at least 3 different components that can 
be distinguished in the image reconstruction process: (f) modulated 
flux F moc i(t) that is modulated by the rotating grids and originates 
from resolved flare sources, (2) unmodulated flux F unm (t), which shows 
no grid modulation but varies proportionally to the smoothed flare 
flux, probably produced by scattered flare photons or by overresolved 
sources that are larger than the pitch of the coarsest grid, and (3) 
non-flare background flux i^ a cfc(0 that may originate from cosmic X- 
ray background, magnetospheric particles, Earth albedo, spacecraft 
radioactivity, or electronic readout noise (Fig. f). Background modeling 
can be performed with the Pixon and Forward-Fit algorithms. We find 
an (unexplained) unmodulated flux in the order of ~ f 5%, in excess of 
the preflare background rate (« 1%). Although the thin attenuator was 
inserted during this flare, it has probably no effect on the unmodulated 
flux because it is located after the grids. 
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Figure 2. Image comparisons of Forward-fit, Pixon, Clean, MEM, and MEM- Vis 
images reconstructed in 4 energy bins. The contour levels are at 10%, 20%, 90% 
of the maximum flux. The location of the northern flux maximum is indicated with 
a dotted cross- hair. The dotted circle has a radius of 18", which is used for halo/core 
flux ratio measurements. 



2.4. Positional Comparisons (§A4 on CD-ROM) 



Images reconstructed in 4 energy bands with all 5 algorithms are shown 
in Fig. 2. The source morphology is a simple double footpoint geometry, 
which is consistently reproduced by all 5 algorithms. MEM- Vis in its 
current implementation does not take corrections for time variability 
into account, and thus seems to have a compromised image quality 
for this case (integrated over a time interval of 40 s). We measure the 
centroid location (xn,Un) of the northern flare source with sub-pixel 
accuracy by parabolic interpolation at the source peak. Averaging these 
source positions among the different image algorithms we find a stan- 
dard deviation in the order of 0.5 pixels (with pixel size of 1"), which is 
close to the accuracy expected from simulations for the corresponding 
count rates. Only MEM- Vis shows occasionally larger deviations, very 
likely to be a side-effect of the compromised image quality. 
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2.5. Amplitude and Width Comparisons (§A5 on CD-ROM) 

The reconstructed source widths w vary to a larger degree between 
individual imaging algorithms. Compared with the spatial resolution of 
the finest used grid (w res pa 11.8" for grid#4), we find that the width 
ratio of the reconstructed nonthermal sources at 25 keV is w/w res = 1.2 
for Forward-Fit, w/w res = 0.7 for Pixon, w/w res = 1.7 for Clean, 
w/w res = 3.3 for MEM, and w/w res = 0.2 for MEM-Vis. The Pixon 
and Forward-Fit algorithms yield compatible widths, while those from 
MEM are significantly larger, and those from MEM-Vis are system- 
atically narrower (probably a side-effect of compromised image quality 
for time- variable sources). The larger widths obtained from Clean maps 
result from the convolution with the Clean beam. Apparently, MEM 
does not converge to the best solution when there are too many degrees 
of freedom. 

The peak flux amplitude a of an image map is roughly reciprocal to 
the square of the source widths, i.e.m a oc 1/w 2 , because the flux per 
feature (F pa aw 2 ) is well-conserved. Because of this strong sensitivity 
to the reconstructed source width, it is not unusual that the peak flux 
differs up to an order of magnitude between different image algorithms, 
depending on the convergence behavior of the reconstruction algorithm. 
For imaging spectroscopy, therefore, it is important to sum over an 
entire feature (such as a flare footpoint) rather than to examine the 
flux of a source on a pixel by pixel basis. 

2.6. Halo/Core Flux Ratios (§A6 on CD-ROM) 

We measure the flux ratio between a core (C) region (inside a radius 
of 1.5 the resolution of the finest grid) and the outer halo (H) region. 
This H/C flux ratio reflects how many photons are spread in the image 
outside the real sources. It is found to be low for Forward-fit, Pixon, and 
Mem-Vis, in the order of a few percent, while Clean shows a somewhat 
larger ratio of H/C £ 30% due to the convolution with the clean beam 
and the addition of clean residuals. MEM exhibits a large noise ratio 
of H/C $ 160%, a clear indication that the image algorithm does not 
converge properly in some cases, in particular for small pixel sizes (1") 
and large field-of- views, which corresponds to a larger degree of freedom 
in the image reconstruction. 

2.7. Spatially Resolved Spectra (§A7 on CD-ROM) 

We fit a two-component (thermal plus non-thermal powerlaw) spectrum 
to the image fluxes obtained in 8 different energy channels. Compar- 
isons of the best-fit parameters with those obtained with the spec- 
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tral fitting software SPEX show agreement between the non-thermal 
fluxes and power-law slopes within a few percents for all imaging algo- 
rithms. However, we emphasize that spatially-resolved spectra should 
only be extracted from images by integrating over an entire source 
structure (or image feature), never from single pixels (because of the 
amplitude-width uncertainty, see §2.5). 

2.8. Grid Selection (§A8 on CD-ROM) 

The grid selection determines the width of the reconstructed sources in 
the case of unresolved sources. We found that the reconstructed source 
width indeed scales (within a factor of 2) with the spatial resolution 
of the finest used grid for most of the image algorithms (except for 
MEM-Vis, which over-resolves, and for MEM, which under-resolves the 
sources). 

2.9. Pixel-size Comparisons (§A9 on CD-ROM) 

The image reconstruction should provide a solution that is independent 
of the map parameterization, such as the field-of-view or pixel size, 
provided that the sources are resolved. Only in the case of unresolved 
sources, the reproduced source width represents an upper limit specified 
by the pixel size. We reconstruct the images with pixel sizes of 1" , 2" , 
and 4" and find indeed an invariant solution for 4 algorithms, with a 
width that is related to the resolution of the finest used grid. Only 
MEM- VIS shows a different behavior, yielding a source size that scales 
with the pixel size rather than with the resolution of the finest used 
grid. 

2.10. Field-of-View Comparisons (§A10 on CD-ROM) 

The selected field-of-view is found not to affect the quality of the im- 
age reconstruction for Forward-fit, Pixon, and Clean. For MEM and 
MEM-Vis, the quality of the image reconstruction seems to improve for 
minimal field-of-views that closely encompass the sources, apparently 
because this minimizes the number of degrees of freedom in the image 
reconstruction. 

2.11. Time Interval Comparisons (§A11 on CD-ROM) 

Since the time interval is proportional to the number of photons used 
in the image reconstruction, the quality and signal-to-noise ratio is 
expected to improve for longer time intervals. We find that the signal- 
to-noise ratio indeed improves with longer time intervals for Pixon, 
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Clean, and MEM-Vis. Forward-fit is a noise- free model by definition, 
so that a noise-to-signal ratio is not defined. For MEM we find no 
noise-to-signal improvement for longer time intervals, when the number 
of degrees of freedom are too large (for small pixel sizes and large 
field-of- views) . 

3. CONCLUSIONS 

This study demonstrates that most available RHESSI image algorithms 
have a reasonable photometric accuracy for the total flux contained 
in the images, but differ substantially in the geometry of the recon- 
structed sources. The results can be summarized as follows: Three 
algorithms (Forward-Fit, Pixon, Clean) converge to the best solution 
within the constraints of the spatial resolution of the finest used grids, 
within the expected uncertainty of the signal-to-noise dictated by pho- 
ton statistics. The current implementation of the MEM algorithm does 
not converge well for cases with a large number of degrees of freedom 
(for small pixel sizes and large field-of- views). The current version 
of MEM-Vis does not correct for time variability and thus yields a 
compromised image quality for long-duration time intervals. Neverthe- 
less, the user can produce consistent and high-quality images with the 
MEM and MEM-Vis algorithms by exercising a careful choice of image 
parameters. 
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